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ABSTRACT 

We use sample of 813 Lyman-break galaxies (LBGs) with 2.6 < z < 3.4 to perform 
a detailed analysis of the redshift-space (z-space) distortions in their clustering pattern 
and from them derive confidence levels in the [T^,/3(z = 3)] plane. We model the 
z-space distortions in the shape of the correlation function measured in orthogonal 
directions, £(<r, 7r). This modeling requires an accurate description of the real-space 
correlation function to be given as an input. From the projection of £(er, it) in the 
angular direction, w p (a), we derive the best fitting amplitude and slope for the LBG 
real-space correlation function: r = 4.48to'is h~ 1 Mpc and 7 = 1.76l u ; u g = 
(?V r o)~ 7 )- A comparison between the shape of £(s) and w p (a) suggests that £(r) 
deviates from a simple power-law model, with a break at ~ 9ft. _1 Mpc. This model 
is consistent with the observed projected correlation function. However, due to the 
limited size of the fields used, the w p (a) results are limited to <r< 10 /i _1 Mpc. Assuming 
this double power-law model, and by analysing the shape distortions in £(<r, it), we 
find the following constraints: f3(z = 3) = 0.151^°, ty n = 0.35±g;||. Combining these 
results with orthogonal constraints from linear evolution of density perturbations, we 
find that (3{z = 3) = 0.25±g;g|, Q? m = 0.55±°;^. 
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1 INTRODUCTION 

To compute the distances to high-redshift objects, one can 
assume that their observed redshifts are only due to the 
expansion of the Universe. Then, for a given assumed cos- 
mology, a relation between comoving distance and redshift is 
easily found. However, the observed redshifts do not exclu- 
sively depend on the expansion rate of the Universe, as they 
also include contributions from the motions of the galax- 
ies in their local rest-frame, superimposed on the systemic 
flow due to the Hubble expansion. If these contributions are 
not taken into account, the estimated distances are said to 
be measured in redshift-space (z-space), and they do not 
reflect the true, comoving distances to the galaxies. If the 
distances could be estimated considering the effects of pe- 
culiar motions, then these would be said to be measured in 
real-space. 

Peculiar velocities introduce distortions in the measured 
clustering pattern. If a given set of galaxies has a spherically 
symmetric clustering pattern in real-space, but the galaxies 
have a very large velocity dispersion, then the clustering 



signal measured in z-space will be smeared along the line-of- 
sight. These features are commonly referred to as "fingers of 
God". Also, gravitational instabilities at large scales cause 
a coherent infall of the galaxies into the potential well of 
higher density regions. 

Besides the z-space distortions, geometric distortions 
in the clustering pattern also occur if a wrong cosmology is 
assumed to convert the measured redshifts into comoving 
distances. This is because the clustering measured along the 
line-of-sight has a dependence on cosmology that differs from 
that of the clustering measurements in the angular direction. 

The goal of this work is to use dynamical and geometric 
distortions to study the large scale coherent infall and put 
constraints on cosmology, from a sample of ~ 800 Lyman- 
break galaxies. The infall can be quantified by measuring f3, 
defined as (3 = f2^ 6 /6 (b being the linear bias parameter). 
The relation between the clustering measured in real- and z- 
space, i n one dimension, in the linear regime, was derived by 
iKaiserl (|l987T) . He found that the large-scale infall will simply 
produce an amplitude shift of the correlation function, mea- 
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sured in real- and z-space. The use of geometric distortions 
to stu dy cosmology was pioneered bv lAlcock fc Paczvnskil 
( 1979), who demonstrated that they can be used as a pow- 
erful cosmological test for a non-zero A. 

Lyman-break galaxies have z ~ 3 and, due to their high 
sky density (~ 1 LBG per arcmin 2 ) and scientific interest, 
such as the study of galaxy formation/evolution at early 
times in the history of the Universe or the existence/role of 
feedback mechanisms at z ~ 3, they are an excellent class of 
objects to study z-space distortions. They can be selected 
by identifying the Lyman continuum discontinuity which, 
at the galaxy's redshift, will determine their location in the 
[U — B,B — R] plane. 

The structure of this paper is as follows: in Section 2 we 
introduce the data used in our analysis. We then measure the 
two-dimensional z-space correlation function that we will 
use to fit the z-space distortions (Section 3). Integrating 
this along the redshift direction allow us to study the real- 
space clustering, as discussed in Section 4. Since there is 
not enough information in the sky direction to use this as 
an amplitude input in the fitting of the z-space distortions, 
we compute the spherical average of the two-dimensional z- 
space correlation function in Section 5, where we determine 
the amplitude and shape of the z-space correlation function, 
in 1-D. In Section 6 we determine the constraints that can 
be obtained from the fitting of the z-space distortions and a 
brief description of the method is given. The discussion and 
conclusions from this work are presented in Section 7. 



2 THE DATA 

T he LBGs used in th is work are included in the sample 
of lSteidel et~ai\ J20031) . The survey presented in their work 
comprises 17 fields, covering an area of 0.38 deg 2 , and the 
total number of 2.67 < z < 3.25 LBGs spectroscopically 
confirmed is 940. 

The details of the comple t e surv ey and the data used 
can be found in ISteidel et al\ l)2003h . [U,G,1Z]ab imaging 
was obtained at several telescopes and used to select the 
LBGs via the Lyman break technique. The spectroscopic 
follow-up was performed at both the Keck I and Keck II 
telescopes. The size of the largest field is 15.5 arcmin and 
the smallest 3.7 arcmin. These limited sizes compromise 
any clustering analysis that depends on the information on 
the sky direction. Seven of these fields contain a background 
QSO, whose spectrum can be used to probe the spatial 
distribution of Lyman a and metal systems. 

The photometric catalogue from which the LBGs were 
selected includes 2347 photometrically selected candidates. 
These have an apparent IZab magnitude limit of 25.5 and 
satisfy the colour criteria: Gab — TIab ^ 1.2, (U{ ni AB} — 
Gab) > {Gab -TIab) + 1.0. 

The determination of the redshift of these candidates 
was, in many cases, derived from the interstellar absorption 
lines of strong transitions, at 1200 — 1700A rest-frame. In 
some cases, the identification of the LBGs was also possible 
by identifying the Lyman-a line. In order to reduce effects 
due to redshift errors that would influence our clustering 
analysis, especially at small scales, we decided to include in 
our sample only the LBGs with class 1 redshift, as defined 



Field Name 


Dimension (arcmin 2 ) 


Number of LBGs 


Q0000-263 


3.69 X 5.13 


15 


CDFa 


8.80 x 8.91 


34 


CDFb 


9.05 x 9.10 


20 


Q0201+1120 


8.69 x 8.72 


21 


Q0256-000 


8.54 x 8.46 


42 


Q0302-003 


15.59 x 15.71 


40 


B20902+34 


6.36 x 6.57 


30 


Q0933+2854 


8.93 x 9.28 


58 


HDF-N 


8.62 x 8.73 


53 


Westphal 


15.02 x 15.10 


176 


Q1422+2309 


7.28 x 15.51 


109 


3C 324 


6.65 x 6.63 


11 


SSA22a 


8.74 x 8.89 


50 


SSA22b 


8.64 x 8.98 


35 


DSF2237a 


9.08 x 9.08 


39 


DSF2237b 


8.99 x 9.08 


42 


Q2233+1341 


9.25 x 9.25 


38 



Table 1. The dimensions of each of the LBG fields in this survey 
and the number of selected LBGs in each field. The field names 
are the same as adopted bv lSteidel et all l2003Fl . When the field 
contains a bright higher- z QSO, the name of the field is the same 
as that of the QSO. Some of the fields are adjacent (CDFa and 
CDFb; SSA22a and SSA22b; DSF2237a and DSF2237b). 



bv lSteidel et al\ l)2003f) . The choice of considering only class 
1 redshift LBGs and the redshift range 2.6 < z < 3.4 leaves 
us with 813 selected galaxies. Tabled shows the number of 
selected LBGs in ea ch of the fields. The field names are the 
same as adopted bv lSteidel et al\ J2003f) . 

The Lyman-a line and the absorption lines are usually 
separated by a factor of a few hundred kms -1 , a feature 
that is often considered as evidence of powerful outflows 
from LBGs. It is prudent to account for this effect, in order 
to have a more precise estimation of the galaxies' redshifts. 
Following Ade lberger et al\ a2003l) we assume a simple out- 
flow model as an explanation for this, where the interstel- 
lar absorption lines are produced "in front" of the outflow 
and hence are blueshifted relative to the galaxy; whereas the 
Lyman-a line is produced in the opposite side of the outflow, 
"behind" the galaxy in the observer's line of sight. Assum- 
ing this simple picture, the systemic redshifts of the LBGs 
can be determined as follows: 

If no absorption features are easily identified and the 
redshift is determined from the Lyman-a line only, then the 
following correction is applied: 



ZLBG 



ZLya 



(1) 



where zlbg is the "corrected" redshift of the galaxy, ZL ya 
the redshift measured from fitting a Gaussian to the Lyman- 
a line's profile, Vtya the mean velocity of the Lyman-a rel- 
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ative to the gal axy's nebular lin e s and c is the speed of 
light. Following lAdelberger et al\ (120031) . we take VL ya = 
310 km s" 1 . 

Similarly, if the redshift is only estimated from the ab- 
sorption lines, then the correction will be: 



ZLBG 



Zabs 



Vabs 
C 



(2) 



where z a bs is the redshift measured from the centroids 
positions of well-defined absorption lines and v a bs is the 
mean velocity of the interstellar absorption lines relative to 
the nebular lines. As we consider the absorption lines to 
be "blueshifted" relative to the galaxy, v a b s will be neg- 
ative. Following Adelbe rger et all <l2003t) . we take v a t s = 
-150 kms -1 . 

In some cases, when both zj^ya and z a b a are measured, 
we apply the following correction (lAdelberger et all l2003): 



ZLBG = 



ZL va + z ab s -O.lUAv + 230 



(3) 



whe re Av = v Lva — v a bs- 

lAdelberger et all ^200A ) . from a subsample of the LBGs 
used in this work, found a value of < Aw >= 614 ± 
314 km s" 1 for the average separation between VL ya and 
v a t s . The "velocity error" of 314 km s" 1 corresponds to an 
uncertainty of ~ 5.6 /i _1 Mpc (comoving separation), that, 
in terms of z-space distortions, produces a similar effect to 
the small-scale galaxy velocity dispersion. 

Fig. Q shows the redshift (zlbg) distribution of our 
sample (solid red line), after computing the corrections 
discussed. This distribution is similar to t he redshift 
distribution presented by ISteidel et all (l2003fl (see their 
Fig. 10). The dashed line is a 4 th order polynomial fit to 
the redshift distribution. 

To quantify the clustering, we must generate a "ran- 
dom" distribution occupying the same volume as the LBGs 
in the survey. In order not to wrongly interpret any com- 
pleteness issue as a feature in the clustering, the random dis- 
tribution must follow the same completeness as the galaxy 
survey. Therefore, we generated a random set of points in 
each the fields of the survey. We generated 20 randoms for 
each LBG in each field. For each of the fields, the redshift 
distribution of this new ensemble is described by the 4 th or- 
der polynomial fit to the LBGs' distribution. Unfortunately, 
we do not have full information about the sky complete- 
ness of the individual fields. Hence, values for the declina- 
tion and right ascension are randomly generated, within the 
field-of-view of each field. The distribution of the number 
of LBGs as a function of angular distance to the centre of 
the field (or, in some cases, the position of the background 
QSO) is consistent with the non-existence of radial gradi- 
ents. It is possible that normalizing the random distribution 
to the number density of galaxies in each field causes some 
clustering features to be "washed out" . If one or more spe- 
cific fields lie directly in regions with particularly high (or 
low) clustering signal, then by generating random sets in- 
dependently of the space density in other fields will cause 
this structure to become unnoticeable. In order to under- 
stand how significantly our measurements are affected by 
this we re-computed the calculations described in Sections 
3, 4, and 5, but using a random set of points generated for 
all the fields simultaneously. The difference between the two 




Figure 1. The solid line is an histogram showing the red- 
shift distribu t ion o f the 813 LBGs taken from the catalogue of 
ISteidel et all i2003t) . The dashed line is a 4" 1 order polynomial 
fit to the distribution. 



results was within the 1 a confidence level. Similarly, by re- 
peating these calculations using different polynomial fits to 
the LBG N(z) as models for the randoms' redshift distribu- 
tion, we concluded that different N(z) models lead to very 
similar clustering measurements. 



3 THE REDSHIFT-SPACE TWO-POINT 
CORRELATION FUNCTION, £(ct,tt) 

Consider the projection of the z-space separation s, between 
two LBGs, along and across the line-of-sight. The z-space 
correlation function can be measured in these two directions, 
which are given by: 



71" = \S2 — Si | 

a= (s2 + s 1 )6>/2 



(4) 
(5) 



where si and S2 are the distances to two different galaxies, 
measured in z-space, and 9 the angular separation between 
them. 

As mentioned above, to estimate the correlation func- 
tion, a random set of points probing the same volume as 
the LBGs must be generated. This ensemble must have all 
the characteristics as the LBGs, such as the sky and redshift 
distributions, although it can not reproduce their clustering. 

Then, the correl ation function £(cr, tt) can be computed 
using the estimator llLandv fc Szalavlll993h : 



< DD(a, tt) > -2 < DR{a, tt) > + < RR(a, tt) > 
< RR(a,n) > ' 



where < DD(a,n) > is the number of LBG-LBG pairs, 
< RR(a, tv) > the number of random-random pairs and 
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a (h 1 Mpc) 



Figure 2. £(<r, 7r) measured for the sample of 813 LBGs, assuming 
a flat, Q^i = 0-3 cosmology. 



goal in this paper is to use these distortions to derive con- 
straints on f3(z = 3) and and draw conclusions on the 
bias and infall of the galaxies at an early stage in the history 
of the Universe, as well as discuss the improvement that can 
be achieved with larger, future LBG surveys. However, the 
study of the real-space correlation function has, just for it- 
self, an obvious interest, as it gives direct information about 
how galaxies cluster, independently of z-space distortion ef- 
fects. 

A picture of the real-space clustering can be obtained 
considering the clustering measured along the a direction, 
since it will not be affected by the z-space distortions. 
This can be obtained by projecting £(<j, ir) along the a di- 
rection, which will give information about th e real- space 
correlation function, £ (r). Following iPeeblesI (Il980f) and 
iDavis fc Peebles! Jl983l) : 



w p (a) = 2 / £(<t, 7r)d7r, 
'o 



(7) 



w p (a) = 2 



< DR(a, 7r) > the number of LBG-random pairs with sepa- 
rations along and across the line-of-sight given by a and tv, 
respectively. 

Fig.|2]shows £(<7, 7r) measured for our LBG sample. The 
shape of the £(<r, n) contours depends greatly on dynamical 
and geometrical effects, whose effects on the clustering can 
be highly anisotropic. 

The LBGs' peculiar velocities lead to distortions in the 
£(<7, 7r) shape, mainly at small scales. The random peculiar 
motions of the QSOs will cause an elongation of the clus- 
tering signal along the 7r direction. The predominant effect 
at large scales is the coherent infall of the LBGs into the 
potential well of overdense regions, which causes a flatten- 
ing of the £(<t, 7r) contours along the ir direction and some 
elongation along a. 

Geometric distortions also occur if the cosmology as- 
sumed to convert the observed QSO redshifts into distances 
is not the same as the true, underlying cosmology of the Uni- 
verse. The reason is because the cosmology dependence of 
the separations along the redshift direction is not the same 
as the one of the separations measur ed in the perpendicular 
direction jAlcock fc Paczvnskilll979Tl . 

Our result does not reproduce the e xtreme elongation 
along the line-of-sight seen in Fig. 20 of lAdelberger et al\ 
( 2003), measured from a subsample of the LBGs used here, 
indicating that the feature was probably due to noise, arising 
from the small number of LBG pairs at small separations. 
T he effects of the velocit y error of ~ 314kms _1 quoted 
bv lAdelberger et al\ (I2003T) are therefore not evident in this 
plot, possibly due to the cancelling of "finger-of-God" effects 
and infall, but more probably due to the effects of noise. 



4 OBTAINING THE PROJECTED 
CORRELATION FUNCTION 

The z-space correlation function consists of a "distorted" 
measurement of the clustering properties of the LBGs. Our 



r£(r) 



\lv 2 — a 2 



(8) 



Fig- EH shows w p (a)/a, obtained from integrating the 
already shown £(cr, n) along the tt direction. To compute 
the errors, we used the Poisson estimate: A£ = (1 + 
£)V 2 / < DD >. The circles are the measured values in 
the present survey. The diamonds are the values found by 
lAdelberger et ali teOOSt) . using a subsample of the LBG en- 
semble for the current work. For a better comparison with 
their results, we used the same values for the upper limit of 
the £(<r, 7r) integration. Hence, this upper limit is the greater 
of 1000 kms " 1 ^ + z)/H(z) and 7a. The fact that the val- 
ues found bv lAdelberger et al\ j2003Tl are systematically be- 
low our results is probably mainly due to differences in the 
random catalogue generated. However, this discrepancy is 
smaller than the error-bars in both results. 

If £(r) can be approximated by a power-law function 
with the form £(r) = (r/ro) ' , it then follows from equation 
IHlthat: 



Up(cr) = 



'rG)r(^)' 



r 5 



(9) 



where T(x) is the Gamma function computed at x. Hence, 
w p (a)/a will also be given by a power-law with the 
same slope as £(r). The best fitting power-law to the 
measured w p (a)/a is represented on the plot by a solid 
line, and it is parameterized by: ro = 4.48^0 ^ /i~ 1 Mpc 
and 7 = 1.76^q'qq. Our £(r) measurement has a higher 
ampli tude and it is a lso steeper than that found by 
(they found r = 3.96±° H, 



lAdelberger et al\ (1200. 



7 = l^to.H)- iFoucaud etal\ ll2003l) . from a sample of 
1294 LBG candidates from the Canada-France Deep Field 
Survey, measured ro for a fixed value of 7 = 1.8 - which is in 
agreement with the slope measured here. Our amplitude is 
still smaller than their w(8) measurements, which indicate 
ro = 5.9t° .l fe^Mpc. 
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Figure 3. The projected correlation function measured for the 
LBG sample. The circles show the result obtained with the respec- 
tive la error bars. The solid line shows the best fittin g power-law 
and th e open diamonds the result obtained by lAdclbcrgc r et all 
fc003fi . 
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Figure 4. The z-space correlation function measured for the LBG 
sample. The circles show the result obtained with the respec- 
tive ler Poisson error bars. The dashed line shows the best-fitting 
power-law to the data and the solid line the best fitting double 
power-law model. The dotted line is the predicted derived 
from the power-law £ (r) model that best describes the w p (<r) data. 



5 OBTAINING THE REDSHIFT-SPACE 
CORRELATION FUNCTION, $(S) 

In order to fit a model to the z-space distortions in £(er, n), 
the correct amplitude of the correlation function must be 
given as an input, in the fitting procedure. Since the fit is 
sensitive to the distortions in the shape of the z-space cor- 
relation function, its correct amplitude must be given as 
an input. Otherwise, the constraints obtained for and 
/3(z) will be such that their values are those needed to com- 
pensate the input £(s), so that the amplitude of the model 
£(<r, 7r) fits the amplitude of £(cr, n) from the data, rather 
then being a good fit to the distortions in £(<r, n). In princi- 
ple, one could use the best-fitting £(r) power-law as an input 
to the £(cr, 7r) model, by decomposing it in two dimensions 
and adding the distortions. However, due to the limited size 
of the fields used in the survey, the behaviour of £(r) at scales 
larger than ~ 10 /i -1 Mpc is unknown and, even considering 
that the power-law approximation is sufficiently good up to 
~ 10 /i -1 Mpc, deviations from a simple power-law model at 
large scales, where £ (a, tt) is fitted, would cause shifts to the 
best-fitting values of /3(z) and O^. The best way of intro- 
ducing the amplitude of the correlation function correctly in 
the £(er, 7r) model is to input a very good description of the 
z-space correlation function's large-scale shape: £(s). 

f(s) reflects the spherical average of £(cr,7r), since s — 
\/n 2 + a 2 . Fig. 2] shows our £(s) measurements from the 
data. 

To avoid underestimating the errors from contamina- 
tion of correlated pairs of galaxies in the same separation 
bin, when the number of pairs is larger than the total num- 



ber of galaxies the ratio 2/ < DD > in the Poisson error 
estimate is replaced by l /N oa i, N aa i being the t otal number 
of galaxies in the survey JShanks fc Bovleil994D . The circles 
show the measured f(s). The error bars represent the la 
confidence level. The dashed line represents the best fitting 
power-law model to £(s). Considering that the line is param- 
eterized by (s/so)~ 7 , then the best fitting values of so and 7 
are given by s = S-llg^ ft _1 Mpc and 7 = 1.711°;™. Given 
the observed deviations from a simple power-law model on 
the f(s) shape (the reduced xL n of the fit is 4.07, with 16 
degrees of freedom), a double power- law model for £(s) was 
also fitted to the data. It was found that the best fitting 
model corresponded to having the two power-laws joining 
at 9/i _1 Mpc. Then, the amplitude of the power-law prob- 
ing the large scales and the slope of both the power-laws 
were fitted. The amplitude of the innermost power-law was 
fixed in such a way that the "break" in the £(s) shape 
was at 9/i _1 Mpc. The parameters characterizing the two 
power-laws are: so = 5.673 /i _1 Mpc and 7 = 1.30!°,;or> 101 
s < 9/i _1 Mpc and s = 7.5±g;| ft _1 Mpc and 7 = 3.291°;^, 
for s > 9/i _1 Mpc. This function is represented by the solid 
line, in Fig. 0] 

The shape of £(s) suggests significant deviations from 
a simple power-law model. The interpretation of this result 
could include the effects of z-space distortions, affecting 
a power-law £(r), whose form is derived from the w p (cr) 
results, as the small-scale random motions of the QSOs 
lead to a deficit of clustering amplitude measured at small 
scales. To confirm this hypothesis, we derived a £(s) model 
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from the best fitting power-law £(r) to the w p (a) results, 
by adding the distortions parameterized by /3(z = 3) = 0.25 
and < vol > 1/2 = 400 km s" 1 . These are like l y valu es for 
those parameters (e.g ., see iFoucaud et al\ i2003l) and 
lAdelberger et al\ J2003)) . To include the distortions in £(r) 
and derive a prediction for £(s), a £(<r, tt) model was derived 
from the f (r) input form, and then this was integrated in 
annuli to obtain £(s). The £(cr,7r) model an d this method 
are described in detail in Ida Angela et al\ (120051) (Model 
1). The dotted black line in Fig. [I] represents the obtained 
£(s) result. This is still a good description of the data, 
within the errors, as it can be seen not only from the 
graph but also from the value of the reduced xL n , which 
is 1.25, for 16 degrees of freedom. Given that these are 
high-redshift galaxies, the z-space distortions in £(s) caused 
by the small-scale random motions are not very significant, 
and hence the observed flattening of £(s) can hardly be 
explained by z-space distortions only. 

Results from several other galaxy and QSO surveys, 
such as the 2dF Galaxy and QSO Redshift surveys, have 
indicated the possible existence of a shoulder, at ~ 8 — 
12 h~ Mpc, in the re spective correlat i on fu nctions (see 
iHawkins eiaTI d2003Tl . Ida Angela et al\ (12005!) '). These re- 
sults are also backed by the shape of the correlation function 
as suggested from CDM model pr edictions, or Halo Occu - 
pation distribution (HOD) models (Ida Angela et all |2005), 
iTinker et ail J2005l) i. 

The double power-law model fitted to the data is a very 
good representation of the £(s) results (the reduced Xmin °f 
the fit is 0.66). The fact that the single power- law model for 
£(r) also represents a good fit to the data, means that we 
can not prefer the double power-law model solely from the 
results of this work. However, when we take into account 
results from other surveys and theoretical models, the dou- 
ble power-law £(s) model is more likely to be closer to the 
correct shape. Future LBG surveys will probe the projected 
correlation function at large separations. Even if these re- 
sults, just from de-projection of w p (a), do not lead to a con- 
clusive measurements of the real-space correlation function, 
then their combination with large-scale statistically reliable 
measurements of £(s) will make it possible to thoroughly 
study the £(r) shape, strongly motivating future wide field 
spectroscopic surveys of LBGs. 

The double power-law £(s) model can therefore be used 
as the input for the amplitude of the 2-space correlation 
function, provided that the fit is only performed on scales 
where the non-linear distortions caused by the random pe- 
culiar motions have a negligible contribution. As it can be 
seen from the dotted black line, which represents the £(s) 
obtained from a simple power-law £(r) model, this is clearly 
not a problem at scales greater than 6 ft -1 Mpc, where the 
function follows very closely a power-law form (assuming the 
value < w 2 z > 1 ' /2 = 400 km s" 1 ). 



6 CONSTRAINTS ON (3 AND O. M FROM 
REDSHIFT-SPACE DISTORTIONS 

Once a model describing how the amplitude of the correla- 
tion function varies with radial separation, in average, then 



the higher-order distortions observed on £(<r, tt) can be mod- 
eled. 

The first step in modeling the z-space distortions is to 
build a £(<7, 7r) model given an input £(r). The latter can be 
obtained from a £ (s) form (7 = 1.30t° 6 7 , for s < 9 /i _1 Mpc; 
7 = 3.29t ?A, for s > 9/i _1 Mpc), if the £(ct,7t) fitting will 
be performed at scales where the main contribution for the 
distortions comes from the coherent infall, and hence the 
input £(r) can be obtained by: 



1 + 



(10) 



The description of the z- space distortion model used 
here is present in section 7.3 of Ida Angela et all J2005fl . Fit- 
ting this model to the measured f (cr, tt) will allow constraints 
on /3(z — 3) to be drawn, provided that a given value for 
the velocity dispersion is assumed. 

Constraints on the density parameter 0,^ are obtained 
from including geometric distortions in the £(<r, tt) model 
and comparing these with the distortions present in £(cr, tt) 
measured from the data. It is useful to make some definitions 
before continuing describing the cosmology f itting through 
geometric distortions. Following iHovle et aZ| (120021) . let the 
true, underlying cosmology of the Universe be the true 
cosmology, the cosmology used to build the model £(o", 7r) 
the test cosmology, and the cosmology assumed, both in the 
model and the data, to measure the correlation function, 
the assumed cosmc 



In the linear regime, the correlation function in the as- 
sumed cosmology will be the same as the correlation function 
in the test cosmology, given that the separations are scaled 
adequately. For £(cr, 7r): 



£t(°"t, TT t ) = £a(o" a , TT a ), 



(11) 



where the subscripts t and a refer to the test and assumed 
cosmology, respectively. 

The relation between th e separations a and tt in the two 
cosm ol ogies is the f ollowing feallinger. Peacock, fc Heavensl 
Il996|) . llHovlell200Cr> : 



(12) 
(13) 

(14) 
(15) 



f Bt 

0~t = /±Cr a = — <7 a 

, At 

where A and B are defined as follows: 
A C 1 



B 



H ^0+^(1 + 3)3 

dz' 



H0j y/nl + fl^il+Z')* 



if the both true and assumed cosmologies correspond to spa- 
tially flat Universes. 

If the same cosmology is assumed in the data and the 
model, then the observed shapes of the measured and mod- 
eled £(cr, it) will be the same when the test cosmology is 
the same as the true, underlying cosmology of the Universe. 
Therefore, constraints on f3(z = 3) and are obtained 
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from simply doing a \ 2 fit between the data and the a se- 
ries of different £(cr, tt) models, derived with different test 
cosmologies and values of (3(z = 3). 

The steps taken for this fitting procedure are as follows: 

• Assume a cosmology and measure £(s), w p (a), £(<j, tt). 

• Take a model for the z-space correlation function, e.g. a 
double power-law. This model should be a good description 
of the observed data at those scales not significantly affected 
by velocity dispersion and included in the £(<r, tt) fit. 

• Choose a pair of test values of f2^i an d P(z). 

• The model for £(s) is a good description for the data 
in the assumed cosmology. What is actually needed at this 
stage is a £(cr, tt) model in some test cosmology, hence the 
correct input for this model is £(s) in that same test cos- 
mology. Since, in the linear regime, £t( r t) = £ a {r a ), one has 
only to compute the real-space separation in the assumed 
cosmology to get £t(r t ). The relation between the real- and 
z-space correlation functions i n any cosmology is given by 
Eq. QT]1 r t is given by r t = yof + (tt* — Wz/H t )' 2 and the 
relation between r a and rt is: r a = rt/(/^/||) 1//3 . 

• Using that model for £t(r t ), compute £t(o~t, Ttt)- Then, 
include the geometric distortions by scaling £t(at,TTt) back 
to the assumed cosmology, in a similar way as described in 
the previous step. To get £ a {a a , TT a ), one needs to scale the 
separations at and 7Vt to a a and Tv a , using equations If 21 and 

EH 

• Adding the effects of large-scale infall not only intro- 
duces distortions in £(cr, tt) but also shifts the amplitude of 
the correlation function, by an amount that depends on the 
value of /3(z) taken. Since the amplitude of the spherical- 
averaged z-space correlation function is introduced as an 
input to the model, it always remains the same, whatever j3 
and il^i are used as test values. This guarantees that the fit 
is being made to the distortions in £(cr, 7r), since the aver- 
aged amplitude remains the same for whatever combination 
of [3 and fi^. 

• For the best fitting value of this amplitude factor, de- 
termine the x 2 value for the fit of this model to the data. 

• Repeat this procedure for a different combinations of 
D,° m and /3(z). 

The number of degrees of freedom in the \ 2 St is the 
total number of bins where £(<r, tt) from the model is fitted 
to the data minus the number of free parameters. If the fit is 
to fi^ and /3(z), the number of free parameters will be two. 
Although it can be argued that the £(a, tt) bins may not be 
independent, the acc uracy of the se errors is supported by 
iV-body simulations iHovle][2000 |). The velocity dispersion 
was fixed to 400 km s" 1 (lAdelberger et al\\2003\) . 

The result of doing this fit in the present data is shown 
in Fig. 

The shaded regions and the solid line refer to the con- 
fidence levels obtained from fitting the z-space distortions. 
These constraints alone correspond to best fitting values of 
/3(z = 3) = O.I5±^° and fi^ = 0.35t° H- The dashed lines 
are the la and 2a confidence levels obtained from clustering 
evolution, from applying linear growth theory and predicting 
the values of (5(z — 3) for different cosmologies, using as an 
input t he value of at z ~ .1 , obtained from the 2dFGRS 
survey (|li^rwkin^^^Z|j2b03) . This method is described in 
detail in Td£^A^igek^^2HT20o"5F) . The dotted line represents 
the la two parameter joint confidence level. The best fitting 
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Figure 5. The confidence levels on the [f2jL , (3(z)] plane from fit- 
ting the z-space distortions in £(<r, tt) (grey scale and solid lines). 
The dashed lines show the la and 2cr confidence levels obtained 
from linear growth theory, considering the value of /3 for local 
galaxies surveys. The joint la confidence level is given by the 
dotted line. 



values are (3(z = 3) = 0.25l° °g and O.^ = 0.551°;^. Fig. 
|S| shows the £(a, tt) data (dashed line) and the best fitting 
f (cr, tt) model (solid line), obtained from the joint constraints 
on /3(z = 3) and A visual comparison between the solid 
and the dashed contours also shows that the model is a good 
description of the £(<r, tt) data. 



7 DISCUSSION AND CONCLUSIONS 

Here we have used the z-space distortions on the cluster- 
ing pattern of a sample of LBGs to derive constraints on 
0(z — 3) and ft^. The method used requires the spherical 
average amplitude of the correlation function to be accu- 
rately known and given as an input to the model. Due to 
the small size of the fields used in the survey, the cluster- 
ing measured across the line-of-sight is not probed at scales 
larger than I0/i -1 Mpc. For this reason, and given that at 
z — 3 only the smallest scales are significantly affected by 
non-linear distortions caused by the random motions of the 
galaxies, the input amplitude of the spherically averaged 
correlation function should be obtained from £(s). 

The £(s) shape suggests that the real-space correlation 
function of the LBGs might deviate significantly from a sim- 
ple power-law model, and these deviations are not evident 
in the projection w p (a). Future LBG surveys where larger 
scales across the line-of-sight are probed, will allow a coher- 
ent picture of the clustering to be drawn from the w p (a) and 
£(s) measurements, similarl y to what has been do ne for the 
2dF QSO Redshift Survey (Ida Angela et al.ll2005l) . 

Once a suitable model for describing the amplitude of 
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a (h 1 Mpc) 

Figure 6. The best fitting model for £(<t, it), obtained using the 
best fitting values of f3{z = 3) = 0.25±2*06 and Q ™ = °- 55 ±o 16 
(solid line) and £(ct, it) measured from the data (dashed line). 



the correlation function as a function of scale is obtained (a 
two power-law model was found to be a good description in 
the present case), then the z-space and geometric distortions 
can be modeled and constraints on f3(z = 3) and fi^n drawn. 
The combination of these constraints with orthogonal con- 
fidence levels from linear growth of density perturbations, 
indicates that (3(z = 3) = 0.25± ;^ and D.° m = 0.55±£;^. 
These va lues are consistent wi th previous measurements of 
the bias. iFoucaud et al\ (2003) found, for a ACDM cosmol- 
ogy, b = 3.5 ±0.3, using a sample of LBGs from the Canada- 
France Deep Field S urvey. Considering the WMAP results 
JSpereel et al\ ho03). if we take = 0.3, which is within 
our computed error bars, their obtained value of b corre- 
sponds to /3(z = 3) = 0.27, which is well within our derived 
error bars and hence consistent with our results. However, 
results from our f (a, n) fits strongly relies on the amplitude 
and shape of the input correlation function. If our input 
£(s) shape is uncertain, then this may introduce errors in the 
measured constraints of and (3{z = 3). We found that the 
derived constraints are very robust to 1 a changes in the pa- 
rameters describing the double power-law £(s) model. How- 
ever, changes to the assumed £(s) (or £(r)) shape, will lead 
to more significant changes in the derived Q.^ and (3(z = 3) 
constraints. Unfortunately, with the present data it is dif- 
ficult to measure the real-space correlation function. This 
handicap is basically due to two factors: the size of fields 
used; and the number of LBGs in the survey. Even though 
the latter is, in part, balanced by the high space density 
of these galaxies, the former prevents any real-space clus- 
tering measurement at scales > 10/i _1 Mpc, except via the 
£(s) form. This hampers any attempt to build such a co- 
herent picture of LBG clustering. We believe we have done 
the best that is possible with the current data. It is im- 
portant to stress the advantage of the LBGs high spatial 
density for z-space distortion analyses. The constraints in 
and (3(z) from the £(o", 7r) fitting alone are compara- 



bleto those achieved with the 2QZ sample (Ida Angela et all 
120051) . which includes ~ 20000 QSOs. However, the much 
smaller number of objects, in the present work, is counter- 
balanced by their high-spatial density, which dramatically 
increases the clustering signal. Our results not only reflect 
what can be achieved in cosmological and dynamical con- 
straints from z-space distortions, but they also foreshadow 
future work, based on larger LBG surveys. 
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